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ABSTRACT 

We report the discovery of a relic Giant Radio Galaxy (GRG) J021659-044920 at 
redshift z ~ 1.3 that exhibits large-scale extended, nearly co-spatial, radio and X-ray 
emission from radio lobes, but no detection of Active Galactic Nuclei core, jets and 
hotspots. The total angular extent of the GRG at the observed frame 0.325 GHz, using 
Giant Metrewave Radio Telescope observations is found to be ~ 2.4 arcmin, that cor¬ 
responds to a total projected linear size of ~ 1.2 Mpc. The integrated radio spectrum 
between 0.240 and 1.4 GHz shows high spectral curvature (cAtm GHz ~ a o'24o GHz 
> 1.19) with sharp steepening above 0.325 GHz, consistent with relic radio emission 
that is ~ 8 x 10 6 yr old. The radio spectral index map between observed frame 0.325 
and 1.4 GHz for the two lobes varies from 1.4 to 2.5 with the steepening trend from 
outer-end to inner-end, indicating backflow of plasma in the lobes. The extended X-ray 
emission characterized by an absorbed power-law with photon index ~ 1.86 favours 
inverse-Compton scattering of the Cosmic Microwave Background (ICCMB) photons 
as the plausible origin. Using both X-ray and radio fluxes under the assumption of 
ICCMB we estimate the magnetic field in the lobes to be 3.3 pG. The magnetic field 
estimate based on energy equipartition is ~ 3.5 pG. Our work presents a case study 
of a rare example of a GRG caught in dying phase in the distant Universe. 

Key words: galaxies: active - galaxies: individual: J021659-044920 - galaxies: jets - 
radio continuum: galaxies - X-rays: galaxies 


1 INTRODUCTION 

The radio morphology of a radio galaxy typically consists of 
a core, a pair of highly collimated jets and lobes formed by 
the supply of relativistic plasma from Active Galactic Nuclei 
(AGN; Scheuer 1974; Begelman et al. 1984; Bridle & Perley 
1984). A less common population of radio galaxies with total 
linear size larger than ~ 0.7 Mpc is referred as Giant Radio 
Galaxies (GRGs; e.g. Saripalli et al. 2005). After an active 
phase, which lasts for ~10 7 — 10 8 yr, the AGN activity stops 
or falls to very low level such that the outflowing jets are 
no longer sustained. As a result, the core, jets and hotspots 
on the lobes disappear. However, the lobes can be seen for 
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relatively shorter period of time (10 6 — 10' yr) before they 
disappear due to radiative losses (Murgia et al. 2011). Thus, 
a relic radio galaxy resulting from the cessation of AGN 
activity represents a short-lived final phase of radio galaxy 
evolution. The short-lived phase makes them a rare class of 
objects. 

The switch off of fresh particle injection leads to an 
exponential steepening of the radio spectrum in the ra¬ 
dio lobes, resulting in a much steeper spectrum (a > 1.3) 
than that of normal radio galaxies (Komissarov & Gubanov 
1994). The low surface brightness and very steep radio spec¬ 
trum makes relic GRGs ideally suited to be detected and 
studied at low frequencies. However, despite concerted ef¬ 
forts to search for relic radio galaxies, hitherto, only a hand¬ 
ful of such objects have been reported in the literature (e.g. 
Cordey 1987; Venturi et al. 1998; Murgia et al. 2005, 2011; 
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Hurley-Walker et al. 2015). In this paper, we report the dis¬ 
covery of a relic GRG J021659-044920 using deep 0.325 and 
1.4 GHz observations carried out at the Giant Metrewave 
Radio Telescope (GMRT) and the Very Large Array (VLA), 
respectively. 

The reported relic GRG lies in the XMM- Large Scale 
Structure (XMM-LSS) field. The existing X-ray data from 
XMM—Newton show extended X-ray emission nearly co- 
spatial with the radio lobes. A variety of physical pro¬ 
cesses can give rise to extended X-ray emission such as 
thermal emission from shocks, synchrotron radiation and 
inverse-Compton (IC) scattering of the ubiquitous Cosmic 
Microwave Background (CMB) photons. Indeed, IC scat¬ 
tered CMB photons have been detected in the form of 
diffuse, extended X-ray emission in several radio galaxies 
at redshift z > 1 [e.g. 4C41.17 (Scharf et al. 2003), 3C 
294 (Fabian et al. 2003), 6C 0905+39 (Blundell et al. 2006; 
Erlund et al. 2008), 4C 23.56 (Johnson et al. 2007), and 
HDF-130 (Fabian et al. 2009), 3C 469.1 and MRC 2216- 
206 (Laskar et al. 2010)]. However, given the dearth of relic 
GRGs and their X-ray observations, X-ray properties of relic 
GRGs are poorly studied. For example, only 2 out of 11 relic 
radio galaxies in Parma et al. (2007) and three out of five 
relic radio galaxies in Murgia et al. (2011) sample have X- 
ray studies. For the relic GRG J021659-044920, we use radio 
and X-ray observations to study its properties e.g., morphol¬ 
ogy, dominant energy loss processes, energetics and magnetic 
field strength. 

In Section 2, we discuss the multiwavelength data used 
in our study. In Section 3, we present the radio morphology 
and spectral properties of the source. The X-ray emission 
properties are discussed in Section 4. In Section 5, we discuss 
the connection between the radio and X-ray emission. We 
estimate the magnetic held strength and the spectral age 
of the relic lobes in Section 6. We summarize our results in 
Section 7. 

Throughout the paper, spectral index is defined as S v oc 
u~ a . We use the Wilkinson Microwave Anisotropy Probe9 
cosmology with Ho = 69.32 km s -1 Mpc -1 , = 0.29 and 

LI a = 0.71. All the frequencies are in the observed frame 
unless mentioned otherwise. 


2 OBSERVATIONS AND DATA ANALYSIS 

The relic GRG J021659-044920 centred at RA 02 h 16 m 59 s 
and Dec. —04° 49 1 20”.6 (J2000) is located in a smaller 
subheld of the XMM-LSS held, known as the Subaru X- 
ray Deep Field (SXDF; Ueda et al. 2008). Both SXDF and 
XMM-LSS have excellent multiwavelength data coverage in 
X-ray, optical, near-infra red (IR), mid-IR and radio bands 
(Simpson et al. 2006; Vardoulaki et al. 2008; Mauduit et al. 
2012). The details of the data used in our study are given 
below. 

2.1 Radio data 

2.1.1 0.325 GHz data 

We use a GMRT observation at 0.325 GHz to study the 
low-frequency radio properties of the relic GRG J021659- 
044920. We have carried out deep 0.325 GHz GMRT ob¬ 


servations of the XMM-LSS held covering ~ 12 deg 2 over 
16-pointings with 32 MHz bandwidth. The GMRT observa¬ 
tions are taken in semisnapshot mode to optimize the uv- 
coverage. The GMRT map has average rms noise of ~ 150 
p Jy beam -1 and angular resolution 9.35 arcsec x 7.38 arc- 
sec. In Figure 1 (left-hand panel), we present the grey-scale 
image of the relic GRG studied in this paper. The details 
of the 0.325 GHz GMRT observations will be presented in 
Wadadekar et al. (in preparation). The 0.325 GHz GMRT 
observations in combination with the existing multiwave¬ 
length data in the XMM-LSS held, have already been used 
to study the population of ultrasteep spectrum radio sources 
(Singh et al. 2014) and the radio-FIR correlation for distant 
star-forming galaxies (Basu et al. 2015). 

2.1.2 0.240 GHz and 0.610 GHz data 

The XMM-LSS held was observed at 0.24 and 0.61 GHz us¬ 
ing the GMRT by Tasse et al. (2007). The 0.240 and 0.610 
GHz observations cover 18.0 and 12.7 deg 2 with rms noise 
~2.5 and ~0.3 mjy beam -1 and angular resolutions of 14.7 
arcsec and 6.5 arcsec, respectively. The Integrated hux den¬ 
sities of the relic GRG reported in Tasse et al. (2007) are 
169.1+20.2 and 20.2+2.7 mjy at 0.240 and 0.610 GHz, re¬ 
spectively. In a deeper GMRT observation of the SXDF at 
0.610 GHz covering 0.5 deg 2 with map rms noise of ~ 60 
/rjy beam -1 , at a resolution of 6.8 arcsec x 5.4 arcsec by 
Vardoulaki et al. (2008), the total hux density of the source 
is found to be 43.3 mjy. This clearly shows that the low 
surface brightness emission is missed in the shallow obser¬ 
vations reported in Tasse et al. (2007). We therefore use the 
0.610 GHz hux density given in Vardoulaki et al. (2008) for 
our analysis. 

2.1.3 1-4 GHz data 

A deep 1.4 GHz image of the source, observed using the 
VLA, is available from Simpson et al. (2006) who has sur¬ 
veyed the SXDF held covering 0.8 deg 2 . The VLA obser¬ 
vations were made using B- and C-array conhguration and 
render rms noise of ~ 20 /iJy beam -1 , at ~ 5.0 arcsec x 4.0 
arcsec angular resolution. We show the 1.4 GHz contours in 
Figure 1 (right-hand panel). 

2.1.4 Missing flux density of radio emission 

The total hux density from a large angular size source can 
be underestimated due to under sampling of the uv-plane at 
the shortest uv- distance. This could be signihcant at higher 
radio frequencies. An interferometer with the shortest base¬ 
line Dm in can detect all the hux from angular scales less than 
~ 0.6A/.D m in, provided the wri-plane is densely sampled at 
the shortest spacings. Here, A is the observing wavelength. 
The total angular extent of the GRG is ~2.4 arcmin at 0.325 
GHz, and can be well sampled by baselines having itw-length 
> 1 kA. This is not an issue with both the GMRT and VLA 
observations at 0.325 and 1.4 GHz, as the shortest baseline 
starts from ~ 0.15 and ~ 0.17 kA, respectively. Hence, we 
do not expect any missing hux density at 0.325 and 1.4 GHz 
and all the structures are well recovered. Additionally, we 
compared the hux density at 1.4 GHz with NRAO VLA Sky 
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Figure 1. Left-hand panel: grey-scale image of the relic GRG J021659-044920 at 0.325 GHz observed using the GMRT having angular 
resolution of 9.35 arcsec X 7.38 arcsec. The map rms noise is ~ 150 //.ly beam” 1 . Overlaid are the contours at 0.325 GHz starting from 
5cr increasing in steps of \/2 X a. The dashed grey contours are at 3rx and —2a. Middle Panel: grey-scale image of the soft X-ray band 
(0.3-2.0 keV) EPIC PN from the XMM— Newton. Overlaid are the 0.325 GHz contours as in the left-hand panel. Right-hand panel: the 
soft X-ray band image overlaid with the 1.4 GHz contours starting from 5<r increasing in steps of \/2 X a. The 1.4 GHz observation has 
angular resolution of ~ 5.0 arcsec X 4.0 arcsec and rms noise of ~ 20 //Jy beam” 1 . In both middle and right-hand panel, the X-ray 
image is smoothed with a Gaussian of kernel radius 5.0 arcsec. The marker ‘x’ indicates the position of the host galaxy. 


Survey (NVSS). They were found to be comparable indicat¬ 
ing negligible missing flux density. 

Further, we note that the 0.610 GHz flux density mea¬ 
sured by Vardoulaki et al. (2008) used baselines above 1.5 
kA for their analysis. Thus, the flux density of 43.3 mjy at 
0.61 GHz could be underestimated due to missing flux and 
hence is lower than the actual flux density. 


2.2 X-ray data 

We used XMM—Newton archival data with observation ID 
0112372001. The data were taken in prime-full window mode 
using thin filter on 2003 January 7, and consist of 25.63 
ks EPIC PN and 27.35 ks of EPIC MOS data. EPIC com¬ 
prises a set of three X-ray CCD cameras i.e., two MOS cam¬ 
eras (Metal Oxide Semi-conductor CCD arrays) and one 
PN camera (pn CCD array). The EPIC cameras perform 
sensitive imaging observations in the energy range 0.15— 
15 keV with spatial resolution 4.0 arcsec — 6.0 arcsec and 
spectral resolution (E/AE) ~ 20 - 50 over the telescopes 
field of view of 30 arcmin (Striider et al. 2001; Turner et al. 
2001). The sensitivity of EPIC MOS/PN sharply declines 
above 10 keV resulting in inadequate counts above 10 keV 
for these observations. Figure 1 shows spatially smoothed 
EPIC-PN soft band (0.3 — 2.0 keV) X-ray image of the 
relic GRG. Overlaid are 0.325 GHz (middle panel) and 1.4 
GHz (right-hand panel) radio contours. The soft X-ray im¬ 
age is smoothed with a Gaussian of kernel radius 5.0 arcsec. 
X-ray emission from extended diffuse double lobes of the 
relic GRG are clearly seen and are nearly co-spatial with 


0.325 GHz lobe emission. We describe the details of X-ray 
data analysis in Section 4. 


2.3 Optical, near-IR data and host galaxy 
identification 

Simpson et al. (2006) present optical identifications for all 
radio sources detected in the 1.4 GHz imaging of the SXDF. 
The optical identifications of radio sources is based on the 
ultra-deep B, R, i!, 2 Suprime-Cam/Subaru images. The op¬ 
tical identifications were made by overlaying the radio con¬ 
tours on to a true-colour B, R, 2 image. The relatively high 
resolution of the 1.4 GHz radio map allowed unambiguous 
identification with optical sources as there was typically only 
one source within an arcsecond of the peak of the radio emis¬ 
sion, even at extremely faint optical magnitudes. For sources 
composed of multiple radio components without core emis¬ 
sion, such as our case, an identification was sought at the 
centre of the source. The relic GRG without a radio core 
presented in this study has been assigned a reliable optical 
counterpart at optical position RA = 02 h 16 m 59 s .064, Dec. 
= —04° 49' 20 ,, .85 with spectroscopically measured redshift 
■Zspec ~ 1.325 (see table 3 in Simpson et al. 2006). The iden¬ 
tified host galaxy of the relic GRG is a very red galaxy (R 
— 2 = 2.0) and resides almost exactly midway between the 
radio lobes. Figure 2 shows overlay of 0.325 GHz radio con¬ 
tours on the optical three-colour RGB image obtained by 
combining Subaru B band (blue), R band (green) and 2 
band (red) images. Optical magnitudes of the host galaxy 
in B, V, R, i' and 2 bands from Simpson et al. (2006) are 
listed in Table 1. The host galaxy is also detected in near- 
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Table 1. AB magnitudes of the host galaxy in the optical and IR bands. 


Band 

B 

V 

R 

i' J 

J 

H 

K 

3.6/zm 

4.5/zm 

5.8^m 

8/im 

Magnitude 

25.08 

24.62 

23.96 

22.98 22.02 

20.30 

19.23 

18.13 

17.83 

17.91 

18.66 

19.23 


Note. Optical band magnitudes ( B, V, R, i', z!) are from Subaru data (Simpson et al. 2006), near-IR fluxes (J, H, K) are from the 
UKIDSS-UDS catalogue and mid-infrared fluxes (3.6, 4.5, 5.6, 8/xm) were measured from the SpUDS images within a fixed circular 
aperture of 6.0 arcsec diameter. 



17m 2hl6m 

Right Ascension (J2000) 


Figure 2. 0.325 GHz (in red) and 1.4 GHz (in yellow) radio 

contours overlaid on the three-colour optical RGB image. RGB 
three-colour image obtained by combining Subaru B band (blue), 
R band (green) and z! band (red) images. The extreme red colour 
of the host galaxy is apparent in RGB image. The arrow points 
to the host galaxy. Both 0.325 and 1.4 GHz contours start from 
5 ct and increase in steps of \/2 X a. 

IR (J, H and A" bands) and mid-IR (3.6, 4.5, 5.8 and 8 /im) 
bands observed as a part of UK Infrared Deep Sky Survey 
(UKIDSS) Ultra Deep Survey (UDS; Lawrence et al. 2007) 
and Spitzer UKIDSS Ultra Deep Survey (SpUDS; PI: James 
S. Dunlop) 1 , respectively. 

2.4 Photometric redshift of the host galaxy 

The spectroscopic redshift of the extremely red host galaxy 
provided by Simpson et al. (2006) is based on a low S/N 
spectrum from FOCAS/Subaru. In order to check the re¬ 
liability of spectroscopic redshift, we estimated photomet¬ 
ric redshift using 12 band photometric data ( B, V, R, i! 
and 2 ! band from Subaru, J, H, K band from UKIDSS and 
3.6, 4.5, 5.0 and 8.0 /im from SpUDS). We estimated the 
photometric redshift (z p hot) of the host by fitting template 
SEDs to the photometric data using the publicly available 

1 http://irsa.ipac.caltech.edu/data/SPITZER/SpUDS/ 


photometric redshift estimation code EAZY 2 (Brammer et al. 
2008). We obtain z p hot = 1.26)ilo'o3, similar to the spectro¬ 
scopic redshift, z spe c = 1.325, within la. Our 2 p hot is also 
consistent with z p hot in the range 0.98 — 1.174 reported 
by Vardoulaki et al. (2008) using relatively shallow mid-IR 
data from SWIRE. Our estimated z p hot independently con¬ 
firms the z sp ec estimate to be reliable. We use z spec = 1.325 
for our calculations. 


3 RADIO PROPERTIES 
3.1 Radio morphology 

We detect the two radio lobes in both the 0.325 and 1.4 GHz 
images, while the AGN core, jets and/or hotspots remain 
undetected. The non-detection of AGN core, jets and hot¬ 
spots classifies this source to be a relic radio galaxy (e.g., 
Murgia et al. 2011; Saripalli et al. 2012). 

At 0.325 GHz, the total integrated flux density of the 
relic GRG J021659-044920 is 140±6.9 mjy with north¬ 
ern and southern lobes having integrated flux densities of 
79.8±5.2 mjy and 60.2±4.5 mjy, respectively. In Table 2, 
we list the flux densities of the GRG at various radio fre¬ 
quencies. The total angular extent at 0.325 GHz (within 50 - 
contour) is 2.4 arcmin and corresponds to a projected lin¬ 
ear size of ~1.2 Mpc at redshift 1.325. The projected linear 
sizes of northern and southern lobes at 0.325 GHz are ~ 600 
kpc (1.18 arcmin) and ~ 510 kpc (1 arcmin), respectively. 
At 1.4 GHz the source size is relatively smaller with total 
end-to-end projected linear size of 1.1 Mpc (2.26 arcmin). 
The total integrated flux density at 1.4 GHz is 9.60 ± 0.14 
mjy with northern and southern lobes having 5.8 ±0.1 mjy 
and 3.80 ±0.09 mjy, respectively. Using 5cr rms noise as the 
upper limit for the undetected AGN core, we estimate the 
core-to-lobe flux density ratio to be < 1 and < 5 per cent 
at 1.4 and 0.325 GHz, respectively. This low value of core- 
to-lobe flux density ratio is similar to other relic GRGs (see 
e.g., Hurley-Walker et al. 2015). 


3.2 Radio spectrum and spectral index map 

In Figure 3 (left-hand panel), we show the radio spectrum 
for the GRG J021659-044920 between 0.24 and 1.4 GHz. 
We list the radio spectral indices computed across 0.240, 
0.325, 0.610 and 1.4 GHz in Table 2. The solid line shows 
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Table 2. Radio flux densities and spectral indices of relic GRG J021659-044920. Column 1: flux density at 1.4 GHz measured using 
Simpson et al. (2006) image. Column 2 is 0.610 GHz flux density of the GRG taken from Vardoulaki et al. (2008). Column 3 is 0.325 
GHz flux density measured using our 0.325 GHz image. Column 4 is 0.240 GHz flux density of the GRG taken from Tasse et al. (2007). 
Columns 5, 6, 7 and 8 are spectral indices between 1.4 to 0.610 GHz, 1.4 to 0.325 GHz, 0.610 to 0.325 GHz and 0.325 to 0.240 GHz, 
respectively. All flux densities are integrated within 5er contours. 



1.4 GHz 
(mJy) 

(1) 

0.610 GHz 
(mJy) 

(2) 

0.325 GHz 
(mJy) 

(3) 

0.240 GHz 
(mJy) 

(4) 

«,1.4 GHz 
“0.610 GHz 

(5) 

~,1.4 GHz 
“0.325 GHz 

(6) 

^0.610 GHz 
“0.325 GHz 

(7) 

^0.325 GHz 
“0.240 GHz 

(8) 

North lobe 

5.80T0.10 

— 

79.8i5.20 



1.79i0.11 



South lobe 

3.80T0.09 

— 

60.2i4.52 



1.89i0.12 



Total 

9.60T0.14 

>43.3 

140i6.9 

169.1i20.2 

> 1.81 

1.83i0.08 

< 1.86 

0.6i0.5 


the best-fitting spectrum fitted by a power-law with expo¬ 
nential cutoff. We describe the spectral modelling in Sec¬ 
tion 3.3. It is evident that the radio spectrum is convex and 
steepens sharply above 0.325 GHz. The source shows strong 
spectral evolution from 0.240 to 1.4 GHz. The radio spec¬ 
tral index estimated between 0.325 and 1.4 GHz is a ~ 1.8, 
consistent with fading lobes without injection of fresh parti¬ 
cles from AGN. In the literature, radio spectral curvature is 
computed as the difference between high-frequency spectral 
index and low-frequency spectral index (i.e, Ohigh — aiow). 
Qhigh — a low (-7 1-0 has been used as a criterion to selected 
relic (dying) radio galaxies (e.g., Murgia et al. 2011). With 
this criterion our source having radio spectral curvature 
(«a6io GHz - «a240 gh!) > 1-19, can be categorized as a dy¬ 
ing radio galaxy. In the active phase, radio galaxies generally 
show radio spectra represented by a power-law over a wide 
range of frequencies. However, in the absence of fresh parti¬ 
cle injection during the relic phase, the energy spectrum of 
the synchrotron emitting relativistic electrons develop an ex¬ 
ponential cut-off towards higher energies due to synchrotron 
and/or IC losses. This leads to an cut-off in the high fre¬ 
quency radio spectrum. Therefore, the fading lobes in relic 
GRG are expected to show very steep (a > 1.3) and convex 
radio spectra, characteristic of electron population that has 
radiated away much of the energy (Komissarov & Gubanov 
1994). 


In Figure 3 (right-hand panel), we present the spectral 
index map computed between 0.325 and 1.4 GHz using pixels 
above 3 a at both the frequencies. The 0.325 GHz emission 
extends well beyond 1.4 GHz emission in the lobes. Spec¬ 
tral index varies from 1.4 to 2.5 with gradual steepening in 
the lobes from outer extremities to inner regions towards 
the host galaxy. The average spectral index of the northern 
lobe is ~ 1.8, and that of the southern lobe is ~ 1.9 (see 
Table 2). The spectral index is extremely steep relative to 
typical low-z radio galaxies where it is ~ 0.7 (Blundell et al. 
1999). However, the variation of spectral index with distance 
from the nucleus in the lobes is consistent with the backflow 
model of the lobes (Leahy & Williams 1984; Leahy et al. 
1989). According to this model, pressure in the hotspot re¬ 
accelerates the post-shock material back towards the core, 
creating a backflow. Thus, the part of the lobe closer to the 
host galaxy consists of an older population of cosmic ray 
electrons (CREs) that are backflowed and makes the radio 
continuum spectrum steep in those regions. 


3.3 Energy losses 

The relativistic electrons in the radio lobes undergo energy 
losses that modifies the power-law radio spectrum. The dom¬ 
inant energy loss processes are synchrotron, IC and adia¬ 
batic cooling. Synchrotron and IC losses have the effect of 
steepening the spectrum and leads to a cutoff at higher ra¬ 
dio frequencies (> 1 GHz), while adiabatic losses affect the 
normalization of the spectrum but not its curvature. In the 
following we suppose that the radiative losses are the dom¬ 
inant energy loss process for the relativistic electrons. We 
neglect the effect of the expansion losses as this is a dying 
source and perhaps in equilibrium with external pressure. 

Both, synchrotron and IC losses affect the spectrum in 
similar ways, whereby, the spectrum is smoothly steepened 
at higher frequencies. The spectrum is characterized by a 
break frequency, ia>i, below which the spectrum remains a 
power-law with a spectral index identical to the injection 
spectral index (Qi n j). Above Vhr, the form of the steepening 
depends on the mechanism of particle injection. For steady 
continuous injection of relativistic electrons (Cl model), the 
spectrum steepens by 0.5 (Pacholczyk 1970). For single shot 
particle injection, the spectrum falls off as a power-law with 
index 4ai n j/3 + 1 above the break frequency, assuming con¬ 
stant pitch angle between electrons and magnetic field (KP 
model; Kardashev 1962; Pacholczyk 1970). However, con¬ 
sidering rapid isotropization of the pitch angle distribution 
leads to an exponential cut-off above the break frequency 
(JP model; Jaffe & Perola 1973). 

In Figure 4 we show the expected spectral index be¬ 
tween rest-frame 0.76 and 3.25 GHz (corresponding to 0.325 
and 1.4 GHz in the observed frame) as a function of the 
break frequency in the rest frame (tv, rest) f° r the different 
models described above. We assumed an injection spectral 
index a; n j = 0.5, and produced synthetic radio spectrum for 
varying i/br,rest- The shaded region shows the observed range 
of spectral index in the lobes. Clearly, the Cl model cannot 
produce the observed steepness and would require Qi n j > 1, 
which is generally not the case. The fact that the extended 
radio emission is very bright at low radio frequencies and 
extends well beyond high frequency radio emission together 
with the observed spectral index supports the fact that the 
lobes consists of an old population of electrons that were 
injected into the intergalactic medium (IGM) from a now 
discontinued jet activity. 

From Figure 4, we find that the observed steepening of 
the spectral index can be produced by the JP model. As per 
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Figure 3. Left-hand panel: integrated radio continuum spectrum between 0.24 and 1.4 GHz of the relic GRG. The spectrum is fitted 
using the JP model for synchrotron radiative loss. Right-hand panel: spectral index map of the GRG estimated between 0.325 and 1.4 
GHz using 3 a flux density levels. Overlaid are the 0.325 GHz total intensity contours. 



Figure 4. Expected spectral index between rest-frame 0.76 and 
3.25 GHz as a function of the break frequency at rest frame, 
k'br.rest for the various models of particle injection with cq n j = 0.5. 
The shaded region shows the observed range of spectral index in 
the lobes. 


the JP model, the range of spectral index suggests, i'br.rest 
to be in the range 1—2.2 GHz in the rest frame. We there¬ 
fore model the observed radio spectrum using the JP model. 
However, we note that, the JP spectrum with a single tv 
is better represented for localized regions. The total radio 
emission is a sum of several JP spectra having a range of 
Vbr- Depending on the time since first injection, i.e., the to¬ 
tal age of the source (f s ) and the time elapsed since injection 
switch-off (toff), the z/br lies in the range tv, low and tv,high, 
respectively. Here, f s = t 0 g + t on , where, t on is the duration 
of the active phase. The two break frequencies are related 
as, I'br.high = ^br,iow(4/toff) 2 (see Murgia et al. (2011)). In 
the scenario, where the dying phase is longer than that of 
the active phase, i.e. toff >> t 0 n , the two break frequencies 
overlap and the global spectrum can be represented by a JP 
model consisting of a single break frequency. In our case, the 
spectrum is not sampled enough to determine the two break 
frequencies. We therefore model the spectrum by a single Vhi 
and approximate the JP spectrum as 52 = J 4z/ _ “ inj e _1 4‘' br . 
Here, A is the normalization at 1 GHz and «i n j is the injec¬ 
tion spectral index. Under this assumption, the zq, r is deter¬ 
mined by t 3 . We assume, a typical value of Qi n j to be 0.5. We 
fit the above form for the JP model to the observed radio 
continuum spectra between 0.24 and 1.4 GHz. For our data, 
we find Vbr = 0.57 ± 0.02 GHz, i.e, zzbr.rest = 1-32 GHz and 
is consistent with what is expected from Figure 4. The best¬ 
fitting radio spectrum is shown as the solid line in Figure 3 
(left-hand panel). 
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Energy (keV) 

Figure 5. 0.3 — 10 keV XMM— Newton PN spectrum best fitted 
with an absorbed power-law, where absorbing column density is 
fixed to the Galactic value. Solid red line and crosses (“+’) repre¬ 
sent fitted model and binned data points, respectively. Residuals 
are shown in the bottom panel. 

4 X-RAY PROPERTIES 

4.1 X-ray data reduction 

We performed X-ray data reduction using Science Analysis 
System (sas) version 13.5 and followed standard procedure 
outlined in the XMM—Newton ABC Guide 3 . We processed 
raw Observations Data Files by using SAS tasks emproc and 
epproc, for the MOS and PN data. We use latest calibra¬ 
tion files to obtain calibrated and concatenated event lists. 
These files were filtered for flaring particle background time 
intervals using standard recommended PATTERN, FLAG 
and energy filters to generate clean event lists. From clean 
calibrated data files we extracted source and background 
spectra. Source spectrum was extracted by choosing two el¬ 
liptical regions covering X-ray emission from each of the two 
lobes. The counts are not sufficient to obtain robust spectral 
fits for individual lobes, and therefore, we opted to extract 
source spectrum covering both the lobes. The background 
spectrum is extracted using a source-free annular region in 
the same chip. The response matrix file and the ancillary re¬ 
sponse file were created using the sas tasks R.MFGEN and 
ARFGEN, respectively. The spectra were grouped to a mini¬ 
mum of 20 counts per bin in order to use y 2 fitting statistics. 
The MOS and PN data are consistent with each other but 
PN data are slightly better at higher energies > 2.0 keV and 
therefore, we present spectral fitting using the PN data only. 

4.2 X-ray spectral fitting and analysis 

To study the nature of X-ray emission from lobes we fitted 
the 0.3 — 10 keV EPIC PN spectrum using physically mo¬ 
tivated models. We use the X-ray spectral fitting package 
XSPEC version 12.8.2 (Arnaud 1996) to fit the X-ray spec¬ 
trum. Since the relic GRG has no AGN core emission, we fit 
the 0.3 — 10 keV spectrum with a single absorbed power- 
law of the form, M(E) = Ke~ N ' a ' T< ~ E) E~ r . Here, K is the 

3 https://heasarc.gsfc.nasa.gov/docs/xmm/abc/ 


normalization of the emission at 1 keV, Ah is the absorb¬ 
ing column density between source and observer, er(E) is 
the photoelectric absorption cross-section (ignoring Thom¬ 
son scattering) and F is the photon index of the power-law. 
We get A h = 2.421®; 33 x 10 21 cm" 2 , F = 2.99t 3 ; 33 and 
y 2 /dof = 7.7/10 ~ 0.77. Keeping all the parameters free 
does not constrain T. However, by fixing the absorbing col¬ 
umn density equal to the Galactic value (Ah = 2.52 x 10 2 ° 
cnY 2 ) statistically improves the fit with y 2 /dof = 9.7/11 ~ 
0.88 and better constrains F to be 1.86/lo (41 ■ We use this 
value of F in our further analysis. The 0.3 - 10 keV EPIC 
PN spectrum fitted with absorbed power-law model for fixed 
Ah is shown in Figure 5 and the best-fitting parameters are 
listed in Table 3. 

To investigate the possibility of thermal origin of X-ray 
emission from hot diffuse gas present in the cluster medium, 
we also fitted the 0.3 — 10 keV spectrum with a thermal 
plasma model (mekal). The mekal model represents an 
emission spectrum from hot diffuse gas based on the model 
calculations of Mewe et al. (1986). The mekal model with 
multiplicative Galactic absorption component does not ren¬ 
der a statistically good fit if all the parameters are kept free. 
We tried fixing the metal abundance of the hot gas in the 
cluster medium in the MEKAL model from 0.1 to 1.5 of Solar 
values in steps of 0.1 but found unphysical or unconstrained 
plasma temperatures (see Table 3). Therefore, it is unlikely 
that X-ray emission from the lobes are of thermal origin. In 
fact, the extended X-ray morphology nearly co-spatial with 
the radio lobes precludes it from being due to hot cluster 
gas. Murgia et al. (2012) reported that relic radio galaxies 
residing in cluster environments display large scale diffuse X- 
ray emission having spherical morphology. And, the X-ray 
emission is characterized by thermal emission from hot gas 
present in the cluster. We note that Finoguenov et al. (2010) 
present a list of X-ray emitting clusters in the SXDF but the 
X-ray emission in our source is reported to be mainly from 
radio lobes. Therefore, we conclude that the observed X-ray 
emission is non-thermal in origin and perhaps originates due 
to IC scattering of the CMB photons (ICCMB). 


5 CONNECTION BETWEEN RADIO AND 
X-RAY EMISSION 

The co-spatial nature of the radio synchrotron and the soft 
X-rays emission from the lobes suggests a fundamental con¬ 
nection between the emitting particle populations within the 
lobes. A possible scenario is IC scattering of CMB photons 
by CREs having Lorentz factors, 7 e > 10 3 , wherein the CMB 
photons are upscattered to X-ray wavebands. Depending on 
the magnetic field strengths, CREs with ~ 10 2 10 5 

give rise to synchrotron emission in the radio wavebands be¬ 
tween ~ 10 MHz and 100 GHz. However, due to radiative 
losses, the lifetime of relativistic electrons are inversely pro¬ 
portional to 7 e . Thus, the radiative lifetimes of (highly ener¬ 
getic) radio synchrotron-emitting electrons are shorter than 
the (less energetic) electrons which give IC emission in the 
X-ray band. As a consequence, IC emission always traces an 
older population of particles which may be more diffuse and 
spatially non-coincident with the radio emission. This gives 
rise to smaller extent and larger off-set between the radio 
emission from the lobes at 1.4 GHz and the X-ray emission 
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Table 3. Best-fitting parameters for the power-law and MEKAL 
models. Errors represent 90 per cent confidence interval on all 
parameters. 


Model 

Parameter 

Value 

power-law 

K-norm 

5.29j; 4 - 39 x 10“ 6 


TVh (cm -2 ) (fixed) 

2.52 x 10 2 ° 


r 

i-86tg;S 


X 2 /dof 

9.7/11 


0.3 — 10 keV flux 1 

3.25 x 10~ 14 


2.0 — 10 keV flux 2 

1.50 x 10~ 14 

MEKAL 

Knorm 

7.90+( S 9 x 10“ 5 


kT (keV) 

4 2+15-3 
^• z -2.3 


nH (cm” 3 ) 

1.17 X KT 2 


Abundance 

0.6 


X 2 /dof 

9.27/10 


1 Observed frame unabsorbed flux in erg cm 2 s 1 . 

2 Rest frame unabsorbed flux in erg cm -2 s —1 . 

(see Figure 1). This is because, at 1.4 GHz the emission 
is produced by a population of relatively short-lived higher 
energy electrons as compared to the 0.325 GHz emission. 

The equivalent magnetic field for IC scattering with 
CMB, Sic, is given by ~ 3.25(1 + z) 2 /xG. For the GRG, 
Bic ~ 17.5 /zG is significantly stronger than that of the 
estimated magnetic field of ~ 3.3 pG (see Section 6.1) indi¬ 
cating IC losses significantly dominates over the synchrotron 
losses. This is because the IC cooling time-scale for electrons 
radiating at 1.4 GHz (i.e., 3.25 GHz in the rest frame) in a 
magnetic field strength of 3.3 /xG is ~ 4.8 x 10 6 yrs, signifi¬ 
cantly lower than that of the synchrotron cooling timescale 
of ~ 1.2 x 10 8 yr. Thus, the radio emission from the lobes 
of the GRG are strongly affected by IC losses. 

The power-law nature of the X-ray emission from the 
lobes can arise due to synchrotron emission from highly en¬ 
ergetic electrons having 7 e ~ 2.7 x 10 s at 2 keV (4.65 keV in 
the rest frame). Such electrons have radiative (synchrotron 
+ IC) lifetime of ~ 2.5 x 10 2 yr in a magnetic field of ~ 

3.3 /xG at z = 1.325. This is about four orders of magni¬ 
tude lower than the estimated spectral age of ~ 8 x 10 6 yr 
in the lobes (see Section 6.3). We therefore conclude that, 
the X-ray emission is a direct consequence of IC scattering 
of the CMB photons with the CREs emitting in the radio 
frequencies. 

Under ICCMB, it is expected that the non-thermal X- 
ray spectrum will have similar spectral shape as of the radio 
spectrum if the IC X-ray and synchrotron radio emission 
are produced due to the same population of electrons. This 
implies that the spectral index in the radio should be the 
same as the slope of the X-ray spectrum. In our observations, 
we find that the radio spectral index, a to be in the range 

1.4 — 2.5 (see Figure 3). The X-ray emission photon index 
(r) is related to the X-ray spectral index (axray) as axray = 
F — 1. Our estimated value of T = 1.86tg4 9 corresponds to 
ax ray = 0.861 q 2 2, significantly flatter than that of the radio 
spectral index estimated between the observed frequencies 
(see Section 3.2). 

In terms of energy, electrons emitting between 0.76 and 


3.25 GHz in the rest frame (i.e., 0.325 and 1.4 GHz in the 
observed frame), correspond to the energy range ~ 3.6 — 7.6 
GeV for magnetic field strengths of 3.3 fiG. Thus, 7 e of the 
electrons are ~ 7 x 10 3 and ~ 1.5 x 10 4 , respectively. For 
CMB photons at z = 1.325, having temperature 6.34 K, the 
Planck function peaks at Uhg ~ 6.6 x 10 11 Hz. These photons 
are up-scattered by the CREs in the radio to X-ray frequen¬ 
cies with average frequencies given by (v) « (4/3)7 2 ;xbg. 
Thus, the CREs emitting at rest frame 0.76 and 3.25 GHz 
would up-scatter these CMB photons to ~ 4.3 x 10 19 and 
~2x 10 2 ° Hz, i.e., in the energy range ~175 and 830 keV. 
Therefore, our X-ray observations between ~ 0.3 and 10 
keV, corresponding to ~ 0.7 and 23.25 keV in the rest 
frame, do not probe this steep part of the spectrum. In 
our observed range of X-ray energies, the emission arises 
due to IC scattering from CREs at much lower energies 
(7e ~ 5 x 10 2 — 1.8 x 10 3 ) corresponding to lower radio 
frequencies where the spectral index is expected to be sig¬ 
nificantly flatter than that between rest frame 0.76 and 3.25 
GHz. In this way, the different spectral indices as measured 
in the X-ray and radio observations can be explained. 

6 MAGNETIC FIELD STRENGTH AND AGE 
OF THE RELIC LOBES 

6.1 Magnetic field strengths 

Together with radio observations, the ICCMB X-ray fluxes 
can be used to constrain the magnetic field strength in the 
lobes of radio galaxies (e.g., Erlund et al. 2006). We use the 
formula given in Tucker (1977) to estimate the magnetic field 
in the lobes as, 



Here, F c and F a are the flux of the X-ray and radio emis¬ 
sion, respectively, Tcmb is the temperature of the CMB 
at the redshift of the source, v c is the frequency of IC 
X-ray emission, v a is the frequency of radio synchrotron 
emission and the constants a(n) and b(n) are taken from 
Ginzburg & Syrovatskii (1965) and Tucker (1977), respec¬ 
tively. 

This analysis requires the observed X-ray and radio 
fluxes to be co-spatial. To get the co-spatial flux, we ex¬ 
tracted the source spectrum in a region matching the 5cr 
contours of the 0.325 GHz image, and fitted the spectrum 
with an absorbed power-law model with Ah set to Galactic 
absorption. The best-fitting photon index is found to be T = 
1.77^0 42 corresponding to spectral index of 0.77+Q' 2 g. The 
co-spatial, unabsorbed, 2 - 10 keV rest frame X-ray flux 

is 1.44 x 10 -14 erg s~ 4 cm -2 that corresponds to 0.75 njy. 
For a synchrotron flux density F s = 140 mjy at v s = 0.325 
GHz, a = 0.77 and hv c « 3 keV, we estimate the B to be 
~ 3.3 /zG. 

6.2 Energy equipartition 

We also compute the magnetic field strength assuming 
equipartition of energy between CREs and magnetic field 
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(see Moffet 1975, equation 2.77). We use the 0.325 GHz emis¬ 
sion, where the radio-emitting lobes of the GRG are assumed 
to be cylindrical in shape, ~ 500 kpc long and ~ 180 kpc 
in diameter each. The synchrotron cutoff frequencies are as¬ 
sumed to be z-Wn = 10 MHz and z/ max = 10 GHz. We neglect 
the contribution of relativistic protons to the total energy of 
relativistic particles in our calculations. Note that, the 1.4 
GHz flux density is heavily affected by ICCMB losses. This 
analysis is valid only for synchrotron losses. Therefore, the 
observed spectral index of 1.83 between 0.325 and 1.4 GHz 
would give rise to magnetic field strengths that are overesti¬ 
mated and would also be sensitive to the choice of lower fre¬ 
quency cut-off. We therefore assume a typical spectral index 
for synchrotron losses as 1.0. Using these values, we estimate 
the equipartition magnetic field strength, B eq ~ 3.5 /rG in 
the lobes. 

The equipartition magnetic field strength and the mag¬ 
netic field strength derived using the X-ray and radio emis¬ 
sion agrees well with each other indicating the energy 
equipartition to be valid in the lobes of the GRG. 

The energy equipartition can be independently tested 
by comparing the total energies in the estimated mag¬ 
netic field strength and the relativistic electrons (£ e ) (see 
e.g., Isobe et al. 2011). We find the unabsorbed 2 — 10 
keV rest frame X-ray luminosity of the lobes to be 1.64 x 
10 44 erg s _1 , consistent with typical luminosities of other 
GRGs (Erlund et al. ( 6 C 0905+39; 2008), Laskar et al. (3C 
469.1; 2010)). Using the X-ray luminosity, we estimate the 
minimum energy in the relativistic electrons (£ e ) following 
equation 4 in Erlund et al. (2006), 


£ e 


3 T 44 64 

4 7 e(l +2) 4 


erg. 


( 2 ) 


Here, L 44 is the X-ray luminosity in units of 10 44 erg s ~ 4 and 
7 e is the typical Lorentz factor of the electrons responsible 
for ICCMB. Assuming 7e = 10 3 , we find £ e = 4.2 x 10 59 
erg. The total magnetic field energy (£b) in the lobes for 
B ~ 3.3 fiG, within an cylindrical volume of ~ 180 kpc di¬ 
ameter and a total length of ~ 1 Mpc is ~ 6.4 x 10 59 erg. 
The two energies match within ~ 30 per cent, indicating 
energy equipartition to be valid. Similar equipartition of en¬ 
ergy has been observed for the GRG DA 240 by Isobe et al. 
(2011). Further, we note that the minimum energy in the 
relativistic electrons for the GRG is similar to other power¬ 
ful high redshift Fanaroff—Riley type II radio galaxies (see 
e.g., Erlund et al. 2006). It shows that the ICCMB electrons 
that are longer lived than their radio emitting synchrotron 
counterparts can deposit a significant amount of energy into 
their surrounding IGM. 


of ~ 3.3 fi G and tv, rest = 1-32 GHz, using Equation 3, 
we find the age of the relativistic electrons in the lobes and 
hence the age of the relic to be t s ~ 8 x 10 6 yr, which is 
typical for relic GRGs (Murgia et al. 2011). 


7 SUMMARY 

We report the discovery of a relic GRG J021659-044920 at 
redshift ( z ) ~ 1.325 with extended diffuse X-ray emission, 
nearly co-spatial with the radio lobes. The salient properties 
of this GRGs are as below. 

(i) The radio lobes are best detected at low radio fre¬ 
quencies observed using the GMRT at 0.325 GHz. The total 
angular extent at 0.325 GHz is 2.4 arcmin that corresponds 
to a projected linear size of ~ 1.2 Mpc at the redshift of the 
source. 

(ii) The host galaxy is identified in deep optical (Subaru), 
near-IR (UKIDSS) and mid-IR (SpUDS). It is a red (R - z 
= 2.0) and dusty galaxy that brightens in mid-IR bands. 

(iii) The relic nature of the radio galaxy is evident as the 

AGN core, jets and/or hotspots remain undetected and the 
lobes exhibits a very steep radio spectral index, gh* 

~ 1.4 — 2.5. The 0.24—1.4 GHz radio spectrum of the lobe 
emission is convex and steepens sharply above 0.325 GHz 
due to radiative losses. The spectral index of the lobes, es¬ 
timated between 0.325 and 1.4 GHz, varies from 1.4 in the 
outer edges to 2.5 in the inner regions, suggesting backflow 
of plasma. 

(iv) The X-ray spectrum between 0.3 and 10 keV is best 
fitted by an absorbed power-law model with photon index, 
r = 1 . 86 l 9 ’ 4 f f° r a fixed absorbing column density similar 
to the Galactic value. The comparison of radio and X-ray 
spectral and morphological properties suggests that X-ray 
emission is likely due to inverse-Comptonization of CMB 
photons by low energy electrons compared to electrons ra¬ 
diating at 0.325 — 1.4 GHz frequencies. 

(v) Using the ICCMB X-ray emission, we estimate the 
lower limit for the total energy in relativistic electrons to be 
~ 4.2 x 10 59 erg (for 7e ~ 10 3 ), implying significant feedback 
from GRG into the surrounding IGM. 

(vi) The magnetic field strength estimated using X-ray 
and radio emission and by energy equipartition yield con¬ 
sistent field strengths of ~ 3.5 fiG. This suggests that the 
energy equipartition between magnetic field and relativistic 
electrons to be valid in the GRG. 
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( 3 ) 


Thus, for our estimated value of magnetic field strength 
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